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Experimen  .al  stress  analysis  techniques  were  used  to  study 
wave  and  fracture  propagation  in  rock  and  rocklike  media.  Trans¬ 
ient  birefringent  specimens  of  CR-39,  Plexiglas  and  Homalite- 
100  were  loaded  with  explosive  charges  and  isochromatic  fringe 
patterns  and  crack  propagation  were  recorded  with  a  multiple 
spark  camera  operating  at  a  rate  of  200,000  frames  per  second. 
Dilatataonal  wave,  shear  wave  and  crack  propagation  velocities 
were  measured.  In  the  case  of  CR-39  crack  velocities  were 
measured  up  to  half  the  shear  wave  velocity,  which  represents 

iheatae°rf t^Cai  velocity.  Marble  and  granite  plates  were 

loaded  explosively  and  isochromatic  fringe  patterns  on  bonded 
photoelastic  coatings  and  moir£  fringe  patterns  were  recorded 

1  nnnennna?  and  Whitley  camera  operating  at  rates  from  250, 00C 
to  1,000,000  frames  per  second.  Dilatational ,  shear  and  Rayleigt 
wave  velocities  were  measured.  The  leading  part  of  the  pulse  is 

an?  ShoWS  ?PPrec^ble  attenuation.  The  trailing  part 
goes  into  tension  causing  widespread  tensile  fracture.  The  ve¬ 
locity  of  propagation  of  this  fracture  zone  in  marble  was  nearly 
the°!feci^1  terminal  velocity.  In  the  case  of  Indeed 
cracks  in  marble  and  granite,  the  velocities  of  crack  propagation 

DD  .r«i".‘\,1473  - - 1 


Security  CUs^i  lic«Uon 


ABSTRACT  (Cont'd. 
DD  Form  1473 


were  appreciably  lower  than  the  theoretical  terminal  values. 
Experimental  results  obtained  were  discussed  and  interpreted 
for  their  relevance  to  the  rapid  excavation  process  in  rock. 
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STUDIES  OF  WAVE  AND  FRACTURE  PROPAGATION 
IN  ROCK  MEDIA 


ABSTRACT 

Experimental  stress  analysis  techniques  were  used  to 
study  wave  and  fracture  propagation  in  rock  and  rocklike  media . 
Transparent  birefringent  specimens  of  CR-39,  Plexiglas  and 
Homalite-100  were  loaded  with  explosive  charges  and  isochromatic 
fringe  patterns  and  crack  propagation  were  recorded  with  a 
multiple  spark  camera  operating  at  a  rate  of  200,000  frames  per 
second.  Dilatational  wave,  shear  wave  and  crack  propagation 
velocities  were  measured.  In  the  case  of  CR-39  crack  velocities 
were  measured  up  to  half  the  shear  wave  velocity,  which  represents 
the  theoretical  limit  velocity.  Marble  and  granite  plates  were 
loaded  explosively  and  isochromatic  fringe  patterns  on  bonded 
photoelastic  coatings  and  moir£  fringe  patterns  were  recorded 
with  a  Beckman  and  Whitley  camera  operating  at  rates  from  250,000 
to  1,000,000  frames  per  second.  Dilatational,  shear  and  Rayleigh 
wave  velocities  were  measured.  The  leading  part  of  the  pulse  is 
compressive  and  shows  appreciable  attenuation.  The  trailing  part 
goes  into  tension  causing  widespread  tensile  fracture.  The 
velocity  of  propagation  of  this  fracture  zone. in  marble  was 
nearly  equal  the  theoretical  terminal  velocity.  In  the  case  of 
induced  cracks  in  marble  and  granite ,  the  velocities  of  crack 
propagation  were  appreciably  lower  than  the  theoretical  terminal 
values.  Experimental  results  obtained  were  discussed  and  inter¬ 
preted  for  their  relevance  to  the  rapid  excavation  process  in 
rock . 


IIT  RESEARCH  INSTITUTE 
iv 


CONTENTS 


Section  page 

1.0  INTRODUCTION  1 

2.0  DYNAMICS  AND  FRACTURE  OF  ROCK  3 

3.0  EXPERIMENTAL  METHODS  8 

3.1  Dynamic  Photoelasticity  8 

3.2  Birefringent  Coating  Method  9 

3.3  Moir4  Method  11 

4.0  WAVE  AND  CRACK  PROPAGATION  IN  ROCK-LIKE  MEDIA  14 

4.1  Experimental  Procedure  14 

4.2  Results  and  Discussion  17 

5,0  WAVE  PROPAGATION  IN  ROCK  MEDIA  31 

5.1  Uniaxial  Wave  Propagation  31 

5.2  Wave  Propagation  in  Marble  Plates  38 

5.2.1  Photoelastic  Coating  Experiments  38 

5.2.2  Moir£  Experiments  54 

5.3  Wave  Propagation  in  Granite  Plates  61 

6.0  FRACTURE  PROPAGATION  IN  ROCK  MEDIA  65 

7.0  SUMMARY  AND  DISCUSSION  72 


REFERENCES 


77 


LIST  OF  ILLUSTRATIONS 


Figure  page 

1  Dimensions  and  Loading  of  Transparent  Photo-  15 

elastic  Specimens.  (a)  Internal  Loading 

( j)  End  Loading 

2  Isochromatic  Fringe  Patterns  in  Specimen  No.  7  18 

(1/8  in.  Thick  CR-39;  50  mg  PETN  Charge;  Camera 
Speed  175,000  Frames  Per  Second) 

3  Isochromatic  Fringe  Patterns  in  Specimen  No.  13  19 

(1/4  in.  Thi?k  CR-39;  100  mg  PETN  Charge;  Camera 
Speed  180,000  Frames  Per  Second) 

4  Wave  Propagation  Velocity  as  a  Function  of  20 

Fringe  Order  in  CR-39  Specimens 

5  Final  Crack  Patterns  in  CR-39  Specimens  No.  4  21 

(1/4  in.  Thick,  100  mg  Charge),  No.  6  (1/4  in. 

Thick,  50  mg  Charge),  and  No,  7  (1/8  in.  Thick, 

50  mg  Charge) 

6  Isochromatic  Fringe  Patterns  in  Specimen  No.  8  24 

(1/4  in.  Thick  Plexiglas;  100  mg  PETN  Charge; 

Camera  Speed  180,000  in. /sec) 

7  Wave  Propagation  Velocity  as  a  Function  of  25 

Fringe  Order  in  Plexiglas  Specimen 

8  Isochromatic  Fringe  Patterns  in  Specimen  No.  10  26 

(1/8  in.  Thick  Homalite-100 ;  50  mg  PETN  Charge; 

Average  Camera  Speed  146,000  Frames  Per  Second; 

Time  Interval  23  p.sec  to  130  ^sec  After  Loading) 

9  Isochromatic  Fringe  Patterns  in  Specimen  No.  11  27 

(1/8  in.  Thick  Homalite-100;  50  mg  PETN  Charge; 

Average  Camera  Speed  170,000  Frames  Per  Second; 

Time  Interval  87  p.sec  to  175  p.sec  After  Loading) 

10  Wave  Propagation  Velocity  as  a  Function  of  28 

Fringe  Order  in  Homalite  100  Specimen 

11  Final  Crack  Patterns  in  Homalite-100  Specimens  29 

No.  10  and  11 


1 1 T  RESEARCH  INSTITUTE 
vi 


LIST  OF  ILLUSTRATIONS  (Cont'd.^ 


Figure 

12  Schematic  of  Instrumentation  for  Uniaxial  Wave 
Propagation  Studies 


13 


14 


15 


16 

17 

18 


19 


20 


21 


Page 

33 


34 


35 


Strain  Gage  Signals  in  Marb?e  Bar  Loaded  by 
aiding  Weight  at  One  End,  (a)  Height  of  Drop 
19  in. ,  (b)  Height  of  Drop  38  in.  (Upper  Trace 
is  First  Gage  4  in.  from  Loaded  End:  Lower  Trace 
is  Second  Gage  4  in.  from  First  Gage) 

pt5a/!r|S  Marble  Bar  Loaded  Explosively  at  One 
End  (Upper  Trace  is  First  Gage  4  in.  from  Loaded 

"iL?wer  Tr*cf  fs  Second  Gage  3.91  in.  From 
First  Gage  and  4  in.  from  Unloaded  End) 

^rani-te  Bar  Loaded  Explosively  at  One 
End  (Upper  Trace  is  First  Gage  and  Lower  Trace 
is  Second  Gage;  the  Distances  Between  Loaded  End, 

First  Gage,  Second  Gage  and  Unloaded  End  are  4  in. 
j./j  m» ,  and  4  in.,  Respectively)  37 

Beckman  and  Whitley  Model  189  Camera  39 

Experimental  Setup  for  Recording  Isochromauic 
Fringe  Patterns  in  Photoelastic  Coating  on  Rock 
Specimens  Loaded  Explosively  43 

Isochromatic  Fringe  Patterns  in  Photoelastic 
Coating  on  Marble  Specimen  Loaded  Explosively  on 
the  Edge.  (Camera  Speed:  500,000  Frames  Per  Second; 
Time  Range:  14-44  p  sec;  Specimen  100372-1)  44 

Isochromatic  Fringe  Patterns  in  Photoelastic 
Coating  on  Marble  Specimen  Loaded  Explosively  on 

Specimen '  1201^2 )  3  5°0’°00  FrameS  Per  Sec°"d;  „ 

Isochromatic  Fringe  Patterns  in  Photoelastic  Coat- 
lng  on  Marble  Specimen  Loaded  Explosively  on  the 
Edge  (Camera  Speed:  500,000  Frames  Per  Second; 

No?6 13 173-1) 18-56  ^  Sec  A£ter  Loading;  Specimen 

Isochromatic  Fringe  Patterns  in  Photoelastic 
Coating  on  Marble  Specimen  Loaded  Explosively  on 
the  Edge . . (Camera  Speed:  500,000  Frames  Per 

13173-1) TimG  Range:  16_5^  ^  Sec;  Specimen  No. 


46 


47 


1 1 T  RESEARCH  INSTITUTE 
vii 


Figure 


LIST  OF  ILLUSTRATIONS  (Cont’d,! 


22 


23 

24 


25 

26 


27 


28 


29 


30 


Isochromatic  Fringe  Patterns  in  Photoelastic 
Coating  on  Marble  Specimen  Loaded  Explosively 
on  the  Edge.  (Camera  Speed:  250,000  Frames  Per 
Second;  Time  Range;  70-98  |i  Sec;  Specimen  No. 
13173-3) 

Wave  Propagation  Velocity  as  a  Function  of 
Fringe  Order  in  Marble  Plates 

Fringe  Order  as  a  Function  of  Location  Along 
Vertical  Axis  with  Time  as  a  Parameter  (Speci¬ 
men  013173-1;  Frame  Numbers  are  Marked) 

Semilogarithmic  Plot  for  Separating  Material 
and  Geometric  Attenuation 

Moir£  Fringe  Patterns  Corresponding  to  Vertical 
Displacements  in  Marble  Specimen  Loaded  Ex¬ 
plosively  on  the  Edge  (Camera  Speed:  1,004,500 
Frames  Per  Second;  Ruling:  1,000  Lines  Per  Inch; 
Specimen  41673-2) 


Moir£  Fringe  Patterns  Corresponding  to  Vertical 
Displacements  in  Marble  Specimen  Loaded  Ex¬ 
plosively  on  the  Edge  (Camera  Speed:  508,000 
Frames  Per  Second;  Ruling:  1,000  Lines  Per  Inch; 
Specimen:  41673-3) 

Vertical  Displacement  Along  Y-Axis  as  a  Function 
of  Depth  with  Time  as  a  Parameter  (Frame  Numbers 
Marked;  Camera  Speed:  1,004,400  Frames/Sec; 
Marble  Specimen  41673-2) 

Vertical  Displacement  Along  Y-Axis  as  a  Function 
of  Depth  with  Time  as  a  Parameter  (Frame  Numbers 
Marked;  Camera  Speed:  508,000  Frames/Sec;  Marble 
Specimen  41673-3) 

Vertical  Strain  as  a  Function  of  Depth  with  Time 
as  a  Parameter  (Frame  Numbers  Marked;  Camera 

??f??:or08’000  Frames/Sec;  Marble  Specimen 
41673-3  * 


Page 


48 

50 

52 

53 


55 


56 


58 


59 


60 


Isochromatic  Fringe  Patterns  in  Photcelastic  Coat¬ 
ing  on  Granite  Specimen  Loaded  Explosively  on  the 
Edge  (Camera  Speed:  500,000  Frames  Per  Second; 

Time  Range:  14-52  \i  Sec;  Specimen  100372-2)  62 


1 1 T  RESEARCH  INSTITUTE 
viii 


Figure 

32 


LIST  OF  ILLUSTRATIONS  (Cont'd.) 


Page 

Isochromatic  Fringe  Patterns  in  Photoelastic 
Coating  on  Granite  Specimen  Loaded  Explosively 
on  the  Edge  (Camera  Speed:  500,000  Frames  Per 
Second;  Time  Range:  22-52  p.  Sec;  Specimen 
030173-2)  63 

33  Wave  Propagation  Velocity  as  a  Function  of  Fringe 
Order  in  Photoelastic  Joating  on  Granite  Specimens  64 

34  Marble  and  Granite  Specimens  After  Loading  66 

35  Isochromatic  Fringe  Patterns  in  Photoelastic 

Coating  on  Marble  Specimen  Loaded  Explosively 
Through  a  Notch  (Camera  Speed:  500,000  Frames 
Per  Second;  Time  Range:  32-70  |i  Sec;  Specimen 
20173-3)  67 

36  Propagation  of  Fracture  Zone  Boundary  in 
Marble  Plate  Leaded  Explosively  Through  a 

Notch  68 

37  Crack  Propagation  in  Marble  Plate  Loaded 
Explosively  Through  Metal  Wedge  in  Notch 
(Frames  Shown  are  at  4  p  Sec  Intervals; 

Specimen  31573-3)  70 

38  Crack  Propagation  in  Granite  Plate  Loaded 
Explosively  Through  Metal  Wedge  in  Notch 
(Frames  Shown  are  at  4  p  Sec  Intervals; 

Specimen  31573-2)  71 


1 1 T  RESEARCH  INSTITUTE 
ix 


STUDIES  OF  WAVE  AND  FRACTURE  PROPAGATION  IN  ROCK  MEDIA 
1.0  INTRODUCTION 

The  increasing  needs  for  urban  and  resource  development 
demand  faster,  safer  and  more  efficient  procedures  for  under¬ 
ground  excavation  in  rock.  The  technology  of  underground 
excavation,  especially  in  hard  rock,  is  not  progressing  in  step 
with  the  needs.  Advancement  of  this  technology  requires  know¬ 
ledge  of  the  geologic  conditions  in  the  region  of  the  excavation, 
i.e.,  the  type  of  rocks  and  soil  present,  and  basic  rock  mechaniro, 
i.e.,  in-situ  rock  properties  and  the  state  of  stress  present. 

The  most  important  process  involved  in  rapid  excavation  is  rock 
disintegration.  The  phenomenon  of  rock  fracture  and  its  relation¬ 
ship  to  the  methods  of  excavation  and  rock  properties  is  not  well 
understood. 

All  methods  of  rapid  excavation  in  hard  rock  must  use 
some  form  of  dynamic  or  impact  loading  on  the  rock,  whether  it 
is  by  means  of  an  explosive  or  a  water  jet.  This  type  of  loading 
produces  stress  waves  which  induce  and  propagate  fracture.  The 
whole  process  of  rock  breakage  then  is  a  complicated  interaction 
of  stress  wave  and  crack  propagation.  The  problem  of  wave  propa¬ 
gation  in  elastic  materials  is  well  understood  but  it  has  not 
been  adequately  studied  in  the  case  of  rock  with  its  intrinsic 
properties  of  inhomogeneity,  anisotropy  and  the  presence  of 
discontinuities.  The  problem  of  fracture  propagation  in  rock 
is  even  more  complex  than  that  of  stress  wave  propagation. 

To  understand  better  the  fracture  processes  in  rock,  it 
is  necessary  to  study  its  response  to  dynamic  loading.  Stress 
wave  propagation  and  fracture  initiation  and  propagation  must 
be  studied  and  correlated.  This  requires  the  measurement  and 
correlation  of  stress  wave  velocities,  i.e.,  dilatat ional,  shear 
and  Rayleigh,  and  the  terminal  crack  propagation  velocity. 
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Analytical  approaches  are  often  limited  as  they  cannot 
fully  account  for  material  inhomogeneity,  anisotropy,  inelasticity 
and  defects.  The  most  promising  approach  then,  is  the  experi¬ 
mental  one.  Techniques  such  as  dynamic  photoelasticity,  photo¬ 
elastic  coatings,  moire  and  strain  gages  are  very  useful  and 
suitable  for  testing  rock  and  simulated  rock  specimens. 

The  objectives  of  this  program  are  to  study  wave  and  crack 
propagation  and  their  interaction  in  rock  and  rock-like  media 
using  the  experimental  techniques  above,  correlate  the  results 
with  available  theory  and  interpret  them  for  their  relevance  to 
the  rapid  excavation  process  in  hard  rock. 
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2.0  DYNAMICS  AND  FRACTURE  OF  ROCK 

Dynamic  behavior  of  rock  is  significant  from  the  point 
of  view  of  rapid  excavation.  Methods  of  rapid  excavation  in 
hard  rock  must  use  some  form  of  dynamic  or  impact  loading  by 
means  of  explosives,  water  jet,  etc.  The  phenomena  of  interest 
under  these  conditions  are  stress  waves,  dynamic  stresses  and 
deformations  and  fracture  initiation  and  propagation. 

The  problem  of  wave  propagation  in  isotropic  elastic 
media  has  been  dealt  with  extensively  by  ai  alytical  and  experi¬ 
mental  methods.  As  described  by  Kolsky2  in  his  monograph  on 
stress  waves,  a  dynamic  loading  in  an  infinite  medium  generates 

two  types  of  waves,  a  dilatational  wave  propagating  with  a 
velocity  of 

r  /  a  +  2U  1 

ci  a) 

and  a  distortional  wave  propagating  with  a  velocity  of 

C2  (2) 

where 

cl>  c2  =  dilatational  and  distortional  wave 
propagation  velocities 

>.»  (x  =  Lame  constants 

p  =  density 

In  addition  to  these  two  basic  types  of  waves  a  third 
type  occurs  along  free  surfaces.  This  wave,  referred  to  as 
Rayleigh  wave,  travels  with  a  vel^ity  slightly  lower  than  that 
of  the  distortional  wave.  The  intensity  of  thesp  waves  decreases 
rapidly  with  depth  below  the  surface. 
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Wave  propagation  in  a  medium  is  defined  by  the  loading 
characteristics  (time  and  space  variation  of  loading),  geometry 
of  the  body  and  material  properties.  Wave  propagation  in  rock 
media  is  a  more  complex  phenomenon  due  to  the  intrinsic  proper¬ 
ties  of  rock,  such  as  anisotropy,  inhomogeneity  and  presence  of 
discontinuities  and  flaws.  Thus,  propagation  velocities  may  be 
different  in  different  directions.  Wave  propagation  may  also 
be  frequency  dependent  due  to  both  inelasticity  of  the  material 
and  to  inhomogeneity.  The  litter  becomes  significant  for  steep 
loading  pulses  when  wavelengths  become  of  the  same  order  of 
magnitude  as  the  grains  and  discontinuities. 

Little  work  has  been  reported  on  wave  propagation  in  rock 

3 

media.  Bienawski  has  measured  dilatational,  distortional  and 
Rayleigh  wave  velocities  in  norite,  where  it  appears  that  these 
velocities  do  not  have  the  same  relative  magnitudes  that  would 
be  predicted  theoretically.  This  points  to  the  need  for  knowing 
the  dynamic  properties  of  the  material.  Ricketts  and  Goldsmith^ 
have  measured  dynamic  properties  in  rock  using  a  Hopkinson  bar 
system. 

The  problem  of  fracture  and  crack  propagation  in  rock  is 
even  more  complex  than  that  of  wave  propagation.  Griffith's 
theory  of  brittle  fracture  states  that  crack  propagation  occurs 
when  very  high  tensile  stresses  are  induced  at  the  tips  of  pre¬ 
existing  cracks  or  flaws  in  the  material. ^  This  theory  was 

6 

further  modified  by  McClintock  and  Walsh  to  account  for  the 
closing  of  cracks  in  a  compressive  stress  field,  and  verified 
experimentally  by  Hoek^  who  used  it  to  explain  failure  around 
openings  in  rock. 

According  to  the  Griffith  criterion  fracture  initiation 
takes  place  when 

a  =  °in.  <3> 
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where 


a  =  applied  uniaxial  tensile  stress 
°in<=  value  of  o  when  crack  propagation  starts 
7  =  specific  surface  energy 
E  =  Young's  modulus 
co  ~  of  preexisting  crack 

An  alternate  form  of  this  criterion  given  by  Irwin®  and 
applied  by  Bienawski®  to  rock  is 


o  =  ^/gE/ttc  ' 

where  G  is  the  energy  released  per  unit  crack  surface  area. 
Fracture  propagation  is  stable  as  long  as  there  is  a  definite 
relationship  between  o  and  c.  When  this  relationship  ceases  to 
exist  and  crack  growth  Ac  is  not  controlled  by  the  stress  incre¬ 
ment  Ao,  crack  growth  becomes  unstable.  This  happens  when  the 

energy  release  rate  reaches  a  critical  value  G  .  Then 

c  * 

0  =  °u  =  fcEhcu‘  (5) 

where  ou  -  value  of  stress  when  fracture  propagation  becomes 
instable 


cu  -  corresponding  crack  half-length 

The  critical  value  Gc  is  a  characteristic  property  of  the 
material. 

When  crack  growth  becomes  unstable,  crack  propagation 
velocity  becomes  important.  Crack  propagation  is  a  dynamic 
phenomenon  and  has  been  treated  as  such  by  Yoffe,9  Roberts  and 
Wells  and  Baker.1  In  these  analyses,  the  problem  is  treated 
either  as  one  of  dynamic  elasticity  or,  in  the  case  of  the 
energy  approach,  account  is  taken  of  the  kinetic  energy  of  the 
material  particles  around  the  propagating  crack.  One  result 
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of  these  analyses  subsequently  verified  by  Wells  and  Post, 
is  that  dynamic  stress  distributions  are  not  appreciably  dif¬ 
ferent  from  static  ones.  Another  significant  conclusion  is 
that  the  crack  velocity  reaches  a  terminal  value 

vt  =  V 2tt / k '  Ve/p  1  (6) 

or,  for  v  =  0.25, 

vt  =  0.38  ViTp1  (7) 

or 

vt  -  <8> 

Once  the  terminal  fracture  velocity  is  reached  crack 
branching  or  bifurcation  occurs,  so  that  increasing  energy  is 
released  through  additional  crack  surfaces.  The  attainment  of 
this  terminal  velocity  coincides  with  ultimate  failure  of  the 
material. 

Experimental  investigations  on  crack  propagation  in  glas 

13 

have  been  pioneered  by  Schardin  who  introduced  the  multiple- 
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spark  camera  (Schardin  camera) .  Wells  and  Post  used  dynamic 
photoelasticity  with  transparent  plastic  models.  Fracture  pro¬ 
pagation  in  rock  has  been  studied  by  Hoek  and  Bienawski^ * 
using  high-speed  photography. 

In  addition  to  fundamental  fracture  mechanics  of  rock 
the  whole  phenomenon  of  rock  blasting  has  been  studied  in  the 
laboratory.  Persson^  used  high-speed  photography  and  plexiglas 
and  glass  models  because  of  their  rock-like  fracture  behavior. 

In  a  series  of  scale  model  blastings,  various  fracture  phenomena 
associated  with  cylindrical  boreholes  were  observed  and  cor¬ 
related  with  derived  stresses  in  the  material.  The  boreholes 
were  charged  with  PETN-powder  to  different  densities,  thereby 
realizing  a  variety  of  detonating  pressures.  The  predicted 


fracture  mechanisms  initiated  by  the  gas  pressure  and  influenced 
by  the  reflecting  shock  waves  were  verified.  Both  two-  and 
ttr-ee -dimensional  studies  were  conducted. 

The  interaction  of  reflected  shock  waves  and  radial 
cracks  moving  out  from  a  single  borehole  has  been  investigated 
by  Field  and  Ladegaard-Pedersen.  A  tensile  wave  passing  the 
material  in  front  of  a  crack  tip  will  contrioute  to  the  stress 
field  with  a  tensile  stress  component  normal  to  the  fracture 
plane  and  the  fracture  velocity  will  increase.  According  to 
Field  those  cracks  travelling  in  directions  forming  an  angle  of 
40-80  degrees  to  the  normal  on  the  free  surface  are  accelerated. 
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3 .0  EXPERIMENTAL  METHODS 
3 . 1  Dynamic  Photoelasticity 

In  dynamic  photoelasticity,  as  in  the  static  case,  a 
model  is  fabricated  out  of  a  birefringent  material  and  subjected 
to  dynamic  loadings  similar  to  those  encountered  in  the  prototype. 
The  model  is  viewed  in  monochromatic  polarized  light  and  the 
resulting  isochromatic  fringe  patterns  are  recorded  by  means  of 
appropriate  electronic  and  photographic  instrumentation. 

The  interpretation  of  fringe  orders  is  the  same  as  in 
the  static  case,  i.e.,  the  fringe  order  is  related  to  the  prin¬ 
cipal  stress  difference  by  the  stress-optic  law 

2nf 

°1  •  °2  *  IT2-  W 

where 

fa  is  the  material  fringe  value 

h  is  the  model  thickness. 

The  value  of  f^  is  determined  by  suitable  calibration 
for  a  dynamic  loading  of  similar  duration  and  rate  as  the  one 
used  in  the  model. 

The  separation  of  stresses,  i.e.,  the  determination  of 
individual  values  of  o^  and  c^j  requires  additional  experimental 
data,  such  as  strains.  The  moire  method  of  strain  analysis  has 
been  used  extensively  for  this  purpose.  In  certain  specific 
cases  of  wave  propagation,  however,  principal  stress  separation 
is  possible  without  any  additional  data.  In  the  case  of  a  pure 
dilatational  wave,  principal  stresses  act  in  the  direction  of 
wave  propagation  and  perpendicular  thereto;  however,  deformations 
occur  only  in  the  direction  of  wave  propagation.  Thus,  the 
displacements  are: 

ur  -  f(r),  u0  =  0 
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(10) 


For  a  point  source,  the  incident  P-^  wave  exhibits  polar 
symmetry,  which  makes  possible  the  determination  of  individual 
strains  (and  stresses)  from  purely  photoelastic  data  as  follows: 


e0 


e 


r 


e0 


1  +  V  n  r 
E  K  V 


(11) 

(12) 


Separation  of  principal  stresses  using  only  photoelastic 
data  is  also  possible  in  the  case  of  a  pure  distortional  wave 
(S  type)  since  for  a  state  of  pure  shear  it  is  known  that 

cx  -  -a 2  (13) 

Thus,  it  is  se:'n  that  separation  of  principal  stresses 
is  easily  accomplished  whenever  there  is  a  distinct  separation 
of  P  and  S  waves . 

Most  birefringent  materials  used  in  model  studies  are 
viscoelastic  in  nature  and,  therefore,  exhibit  rate  (or  time, 
or  frequency)  dependence  in  all  mechanical  and  optical  proper¬ 
ties.  Whenever  material  behavior  is  such  that  rate  effects  cannot 
be  neglected,  a  photoviscoelastic  analysis  of  the  fringe  patterns 
is  necessary.  Methods  and  techniques  for  such  an  analysis  have 
been  developed  and  successfully  applied  to  wave  propagation 
problems.  The  stress-optic  law  in  this  case  takes  more  compli¬ 
cated  forms,  involving  a  time-dependent  stress  fringe  value 
which  must  be  determined  by  proper  dynamic  calibration  tests. 

3 . 2  Birefringent  Coating  Method 

While  the  conventional  method  of  transmission  photo¬ 
elasticity  is  not  feasible  with  opaque  rock  materials,  the 
birefringent  coating  technique  is  ideally  suited  for  determining 
the  surface  strains  in  such  materials.  Knowledge  of  the  strain 
field  and  the  stress-strain  response  for  the  rock  enables  evalua¬ 
tion  of  the  stress  field.  Numerous  studies  have  also  been 

NT  RESEARCH  INSTITUTE 


9 


i. 


conducted  to  further  develop  the  method  and  to  establish  the 
accuracy  of  the  technique  under  various  practical  conditions. 

The  birefringent  coating  technique  consists  of  bonding 
a  thin  sheet  of  photoelastic  material  to  the  structure,  in  this 
case  the  rock.  When  the  rock  structure  is  loaded,  the  surface 
deformations  and  strains  of  the  rock  are  transferred  to  the 
borded  coating  and  induce  a  birefringence  in  the  latter.  Tne 
birefringence  is  observed  and  recorded  by  means  of  a  reflection 
polariscope,  where  the  polarized  light  traverses  the  thickness 
of  the  coating  .d  is  reflected  back  from  a  reflective  layer 
provided  at  the  interface  between  the  coating  and  the  structure. 
Within  the  linear  strain-optic  range  of  the  coating,  the  fringe 
patterns  are  directly  related  to  the  surface  deformations  (and 
strains)  in  the  structure  whether  they  are  elastic  or  not. 

While  strain  gages  provide  point -by-point  results,  photoelastic 
coatings  give  full-field  information.  Moreover,  coatings 
provide  instantaneous  visual  observation  of  the  strain  field 
throughout  the  loading  pulse. 

At  any  generic  point,  the  relative  retardation  fringe 
order  n  of  the  photoelastic  coating  is  related  to  the  difference 
in  the  in-plane  principal  strains  in  the  rock  e^,  E£  by  the 
strain-optic  law: 

el  "  e2  =  =  ~7t  (14) 


where 


=  coating  thickness 
=  wavelength  of  light 
=  strain  sensitivity  coefficient 
=  strain  fringe  value  of  coating  material 
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In  those  cases  where  the  material  response  is  linearly 
elastic,  the  observed  birefringence  in  the  coating  is  related 
to  the  principal  stresses  o2  in  the  structure  (rock)  by 

°1  '  °2  =  “Zt  (1  +'" v)  (15) 

where  E  and  v  are  the  elastic  modulus  and  Poisson's  ratio  of 
the  material. 

The  accuracy  of  the  coating  technique  is  influenced  by 

the  mismatch  in  Poisson's  ratio  between  the  rock  and  coating. 

While  the  Poisson's  ratio  for  the  rock  materials  is  lower  than 

that  of  the  coating,  the  mismatch  effect  can  be  minimized  by 

using  thin  coatings.  However,  thin  coatings  produce  fewer 

fringes  as  can  be  seen  from  Eq .  (14) .  The  mismatch  in  Poisson's 

ratio  becomes  significant  over  a  transition  zone  between  a  free 

boundary  of  a  specimen  and  the  internal  region.  Photoeiastic 

coatings  are  ideal  for  detecting  visually  regions  of  plastic 

deformation,  fracture  initiation  and  propagation.  Photoelastic 

coatings  have  been  applied  successfully  to  the  study  of  strain 

19  20 

distribution  and  fracture  around  cavities  in  rock.  ’ 

3.3  Moire  Method 

The  moire  effect  is  an  optical  phenomenon  observed  when 
two  closely  spaced  arrays  of  lines  are  superimposed  and  viewed 
with  either  transmitted  or  reflected  light.  If  the  two  arrays 
consist  of  opaque  parallel  lines  which  are  not  identical  in 
either  spacing  (pitch)  or  orientation,  then  interference  between 
the  two  arrays  occurs  and  moire  fringes  are  produced.  The  con¬ 
cept  offers  a  reasonably  easy  means  of  full-field  in-plane 
displacement  and  strain  analysis.  Being  strictly  geometric  in 
nature,  the  technique  is  independent  of  the  anisotropy  and  the 
stress-strain  response  of  the  material.  The  in-place  displace¬ 
ment  field  is  measured  directly  and  the  strain  fields  are  derived 
by  differentiation  of  these  displacement  fields. 
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In  actual  practice,  one  array  is  photo-printed,  etched, 
or  bonded  on  the  surface  of  the  specimen.  A  transparent  refer¬ 
ence  array  usually  of  the  same  pitch  is  placed  adjacent  to  the 
specimen  and  aligned  parallel  to  the  specimen  array.  As  the 
specimen  is  loaded,  the  specimen  array  will  deform  and  follow 
the  surface  displacements  induced  by  the  loads.  The  reference 
array,  of  course,  does  not  deform  since  it  is  not  loaded  and  as 
a  consequence  a  moire  interference  pattern  forms  when  light  is 
transmitted  through  the  two  arrays. 

Moire  fringe  patterns  can  most  easily  be  interpreted  by 
relating  them  to  the  displacement  field  since  moir4  fringes  are 
loci  of  points  having  the  same  component  of  displacement  in  a 
direction  perpendicular  to  the  grid  lines  of  the  undistorted 
reference  array.  Once  two  orthogonal  displacement  fields  (u  and 
v)  have  been  established  by  using  line  arrays  parallel  to  the 
x  and  y  axes  of  a  specimen,  the  cartesian  strain  components  can 
be  computed  from  the  derivatives  of  the  displacements.  For 
small  strains  these  relations  take  the  form: 

£X  ?iX 

e  = 

y  ay  (16) 

y  =  Su  av 

'xy  ay  ax 


where 


x,y  -  cartesian  reference  axes 
u»v  =  x  and  y  components  of  the  displacement 
ex,ey  =  normal  components  of  strain 
7xy  =  shear  component  of  strain 
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The  displacement  gradients  ?}u/^x  and  ?>vfoy  are  obtained  from 
the  slopes  of  the  two  displacement  surfaces  in  a  direction 
perpendicular  to  the  lines  of  the  reference  arrays.  The  dis¬ 
placement  gradients  Su/^y  and  ^v/^x  are  obtained  from  the  slopes 
of  the  displacement  surfaces  in  a  direction  parallel  to  the  lines 
of  the  reference  array. 

The  fineness  of  the  line  images  required  to  obtain 
satisfactory  information  from  a  moire  pattern  depends  only  upon 
the  magnitude  of  the  deformations.  However,  these  deformations 
are  in  turn  related  to  the  modulus  of  elasticity  of  the  specimen 
material,  the  geometry  of  the  specimen,  and  the  manner  of  load¬ 
ing.  Line  spacings  for  a  given  material  can  only  be  specified 
over  a  broad  range.  The  pitch  of  the  array  should  be  small 
enough  to  give  sufficient  response  even  in  cases  where  the 
material  is  stressed  below  the  yield  point  (in  the  case  of 
elastic  analyses).  On  the  other  hand,  the  pitch  should  not  be 
so  fine  as  to  produce  fringes  that  can  no  longer  be  clearly 
distinguished.  The  moire  method  is  well  suited  for  the  study 
of  transient  displacement  (hence  strain)  fields  and  fracture 
propagation  in  dynamically  loaded  rock  specimens. 
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4.0  WAVE  AND  CRACK  PROPAGATION  IN  ROCK-LIKE  MEDIA 
4 . 1  Experimental  Procedure 

The  materials  used  in  this  phase  of  the  program  were 
transparent  birefringent  materials  to  allow  the  application  of 
dynamic  photoelastic  techniques.  Three  different  materials  were 
used,  Columbia  Resin  (CR-39),  Polyester  Resin  (Homalite  100), 
and  Plexiglas.  The  specimens  were  10  in.  x  6  in.  plates  1/4  in. 
and  1/8  in.  thick.  Thirteen  specimens  were  tested  by  means  of 
explosive  loading.  Most  specimens  were  loaded  through  an 
internal  1/2  in.  diameter  hole  located  along  the  long  center- 
line  2.5  in.  from  one  end  (Fig.  la).  A  1/2  in.  long  radial 
slot  was  sawed  at  one  point  on  the  hole  boundary  along  the 
horizontal  centerline  to  facilitate  crack  initiation  at  that 
point.  Some  specimens  were  loaded  through  a  slot  at  one  end 
(Fig.  lb). 

Explosive  loading  was  produced  by  detonating  50  mg  and 
100  mg  charges  of  Pentaerythritol  Tetranitrate  (PETN) .  This 
explosive  is  readily  available  in  Primacord  and  is  well  suited 
for  dynamic  loading  of  models.  Typical  pulses  produced  by  PETN 
have  a  rise  time  of  approximately  one  microsecond  and  a  mean 
width  of  six  microseconds.  Furthermore,  PETN  explosions  do  not 
produce  too  much  smoke  or  light,  thus  allowing  photography  of 
fringe  patterns  near  the  explosive  source. 

The  PETN  powder  was  dissolved  in  acetone  and  precipitated 
in  ice  water  to  obtain  a  finer  particle  size.  This  powder  was 
then  pressed  into  specially  made  plexiglas  detonator  caps.  The 
cap  with  the  explosive  was  attached  to  an  exploding  bridge  wire 
header.  The  EBW  header  consisted  of  two  terminal  pins  embedded 
in  a  plastic  cylinder.  The  terminal  pin  spacing  is  bridged  with 
0.0015  in.  diameter  gold  wire.  When  a  high  voltage  (2000v)  is 
rapidly  (in  1  p.  sec)  imposed  across  the  header  terminals  the 
gold  wire  explodes  and  initiates  the  detonation  of  the  fine 
grain  PETN  in  contact. 

I 
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EXPLOSIVE 

CHARGE 


EXPLOSIVE 

CHARGE 


Fig.  1  DIMENSIONS  AND  LOADING  OF  TRANSPARENT  PHOTOELASTIC  SPECIMENS 
(a)  INTERNAL  LOADING  (b)  END  LOADING 


Dynamic  photoelastic  fringe  patterns  were  recorded  by  a 
multiple  spark  camera  (Cranz-Schardin) .  This  camera  system  is 
capable  of  photographically  recording  a  time-sequenced  series 
of  16  dynamic  photoelastic  fringe  patterns  at  rates  that  can  be 
varied  between  185,000  and  700,000  frames  per  second.  The  line 
definition  of  the  photographic  records  is  such  that  an  array  of 
approximately  12  lines  per  inch  moving  at  a  velocity  of  100,000 
in. /sec  can  be  resolved. 

The  light  source  is  provided  by  a  4  by  4  array  of  16 
separate  spark  gaps  discharged  sequentially.  The  exact  time  and 
sequence  of  firing  is  sensed  with  a  photosensitive  device  (Lite 
mike)  and  recorded  with  an  oscilloscope.  The  light  from  each 
spark  gap  passes  through  a  circular  polaroid  element,  the  model 
and  a  second  circular  polaroid  element.  It  is  then  collected 
by  a  9-in.  diameter  field  lens  and  is  focused  on  one  of  16 
camera  lenses.  To  obtain  a  narrow  bandwidth  of  light,  gelatin 
filters  (No.  47B)  having  a  maximum  transmission  between  4,200 
and  4,400  angstroms  were  placed  behind  the  camera  lenses.  The 
16  images  obtained  are  recorded  on  a  large  (11  in.  x  14  in.) 
sheet  of  Kodak  Commercial  Film  No.  4127. 

In  operation  of  the  system,  a  number  of  events  were 
coordinated;  the  detonation  of  the  explosive  charge,  the  trig¬ 
gering  of  the  oscilloscope  and  the  triggering  of  the  spark  gap 
bank.  The  sequence  of  events  was  initiated  by  means  of  a 
firing  circuit,  emitting  simultaneously  three  pulses:  (1)  a 
2  kv  pulse  which  exploded  the  bridge  wire  and  detonated  the 
PETN,  (2)  a  20  v  pulse  which  triggered  the  oscilloscope,  and 
(3)  a  20  v  pulse  which  triggered  the  spark  gap  bank  through  a 
time  delay  generator. 

Four  1/4  in.  thick  CR-39  specimens  were  loaded  with 
100  mg  PETN  charges,  and  two  1/8  in.  thick  CR-39  specimens 
were  loaded  with  50  mg  charges.  Two  1/4  in.  thick  Plexiglas 
plates  were  loaded  with  100  mg  charges  and  two  1/8  in.  thick 
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Homalite-100  plates  were  loaded  with  50  mg  charges.  The  resulting 
fringe  and  crack  patterns  were  recorded  by  the  Cranz-Schardin 
camera  at  a  rate  of  180,000  frames  per  second. 

4 . 2  Results  and  Discussion 

Typical  records  of  dynamic  fringe  patterns  in  1/8  in. 
and  1/4  in.  thick  CR-39  specimens  loaded  with  50  mg  and  100  mg 
charges  of  PETN  are  shown  in  Figs.  2  and  3.  The  propagation 
velocities  of  various  fringe  orders  were  computed  from  such 
records  and  plotted  versus  fringe  order  in  Fig.  4.  The  results 
show  some  spread  due  to  the  variation  of  specimen  thickness  and 
amount  of  explosive  charge.  The  wavefront  velocity  obtained 
by  extrapolation  to  the  zero  fringe  order  varied  between  75,500 
in. /sec  and  80,000  in. /sec.  The  velocity  of  the  shear  waves 
reflected  from  the  free  edges  of  the  specimen  was  also  measured. 

It  varied  between  38,700  in. /sec  and  41,500  in. /sec,  the  higher 
value  obtained  from  the  1/8  in.  thick  specimen  (No.  7).  This 
velocity  can  also  be  calculated  from  the  measured  dilatational 
plate  velocity  and  Poisson's  ratio  as  follows: 


where  cg  =  shear  velocity 

=  dilatational  plate  velocity 
v  =  Poisson's  ratio 

For  an  average  value  of  cL  =  78,000  in. /sec  and  v  =  0.35, 

<^s  =  44,400  in. /sec. 

Excessive  crack  branching  was  observed  in  all  specimens 
due  to  the  high  energy  of  the  explosive  charge.  Final  crack 
patterns  for  CR-39  specimens  No.  4,  6,  and  7  are  shown  in 
Fig.  5.  The  crack  propagation  velocities  varied  more  widely 
than  wave  propagation  velocities.  They  were  17,200,  14,000, 
20,000  and  20,500  in. /sec  for  specimens  No.  4,  6,  7,  and  13. 
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Fig.  2  ISOCHROMATIC  FRINGE  PATTERNS  IN  SPECIMEN  NO.  7  (1/8  IN.  THICK 

CR-39 ;  50mg  PETN  CHARGE;  CAMERA  SPEED  175,000  FRAMES  PER  SECOND) 
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Fig.  3  ISOCHROMATIC  FRINGE  PATTERNS  IN  SPECIMEN  NO.  13  (l/<+  in.  Thick  CR-39; 
100  mg  PETN  Charge;  Camera  Speed  180,000  Frames  Per  Second) 
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Fig.  4  WAVE  PROPAGATION  VELOCITY  AS  A  FUNCTION 
OF  FRINGE  ORDER  IN  CR-39  SPECIMENS 
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Fig.  5  FINAL  CRACK  PATTERNS  IN  CR-39  SPECIMENS  NO.  4  (1/4  in 
THICK,  100  mg  CHARGE),  NO.  6  (1/4  in.  THICK,  30  mg 
CHARGE),  AND  NO.  7  (1/8  in.  THICK,  30  mg  CHARGE) 
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Wave  and  crack  propagation  results  for  CR-39  obtained 

from  four  tests  are  summarized  in  Table  1.  It  is  seen  that  the 

maximum  ratio  of  the  crack  propagation  velocity  jo  the  shear  wave 
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velocity  is  0.50,  equal  to  the  theoretical  prediction.  The 
lowest  ratio  was  found  for  the  1/4  in.  thick  specimen  with  the 
50  mg  explosive  charge. 

Tests  similar  to  those  described  above  were  conducted 
with  Plexiglas  and  Homalite-100  specimens.  A  typical  record 
of  dynamic  fringe  patterns  in  a  Plexiglas  specimen  is  shown  in 
Fig.  6.  The  propagation  velocity  as  a  function  of  fringe  order 
is  plotted  in  Fig.  7.  The  wavefront  velocity,  obtained  by 
extrapolation  to  the  zero  fringe  order,  is  96,000  in. /sec.  The 
shear  wave  velocity  was  obtained  by  three  different  ways.  The 
velocity  of  the  incident  shear  wave  visible  in  the  fringe 
patterns  of  Fig.  6  is  53,700  in. /sec,  that  of  the  shear  wave 
reflected  from  the  free  edges  is  50,000  in. /sec,  and  that  cal¬ 
culated  from  the  measured  dilatational  wave  with  an  assumed 
value  of  0.35  for  Poisson's  ratio  is  54,600  in. /sec.  A  very 
limited  amount  of  cracking  was  observed  because  of  the  more 
pronounced  viscoelasticity  and  ductility  of  the  material. 

Dynamic  fringe  patterns  in  1/8  in.  thick  Homalite-100 
plates  loaded  with  50  mg  PETN  charges  are  shown  in  Figs.  8  and 
9.  Figure  8  shows  the  first  stage  of  the  phenomenon  covering 
the  period  between  23  p,sec  and  130  p,sec  after  loading,  and 
Fig.  9  covers  a  later  stage  between  87  ^sec  and  175  ^sec  after 
loading.  The  propagation  velocity  as  a  function  of  fringe 
order  is  plotted  in  Fig.  10.  The  wave-front  velocity  obtained 
by  extrapolation  is  cL  =  89,000  in. /sec.  The  incident  shear 
wave  velocity  is  49,500  in. /sec,  and  the  one  computed  from  the 
measured  dilatational  velocity  is  50,500  in. /sec.  Cracking 
again  was  limited  and  the  maximum  crack  propagation  velocity 
observed  was  13,500  in. /sec,  far  below  the  limiting  theoretical 
value  of  25,000  in. /sec.  Final  crack  patterns  of  the  two 
Homalite-100  specimens  are  shown  in  Fig.  11.  The  pattern  of 
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specimen  No.  10  shows  excessive  crack  branching  with  both  radial 
and  circumferential  cracking  around  the  loaded  hole.  Spalling 
at  the  corners  and  at  the  edge  farthest  from  the  hole  due  to 
reflected  tensile  waves  is  plainly  visible. 


I 
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5.0  WAVE  PROPAGATION  IN  ROCK  MEDIA 


5.1  Uniaxial  Wave  Propagation 


Vermont  marble  and  Cold  Spring  (Minnesota)  granite  were 
the  two  rock  media  investigated.*  Bar  specimens  of  dimensions 
1/2  in.  x  1/2  in.  x  12  in.  were  used  to  determine 
uniaxial  wave  propagation  and  the  dynamic  tensile  strength  of 
these  materials. 


The  stress-wave  technique  has  been  employed  for  this 
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purpose  to  determine  the  dynamic  properties  of  brittle  materials. 

A  slender,  freely  supported  bar  is  impulsively  loaded  at  one  end 
and  strains  or  displacements  are  measured  at  one  or  two  stations 
along  the  bar.  The  wave  propagation  velocity  is  determined  from 
the  travel  time  of  the  pulse  between  two  stations.  This  uniaxial 
wave  velocity  is  related  to  the  dynamic  modulus  of  the  material 
(E)  and  its  density  (p)  by 


(18) 


The  stress  at  any  point  can  be  determined  from  measure¬ 
ment  of  either  the  axial  strain  or  the  particle  velocity  at  that 
point.  The  axial  stress  is  determined  by  the  following  relation¬ 
ships  : 

°x  pco  ST  "  Eex  <19> 

where 

a  =  axial  stress 

=  particle  velocity  at  point  in  question 
ev  =  axial  strain  at  point  in  question 

A 

JL. 

The  granite  was  supplied  free  of  charge  by  the  Cold  Spring 
Granite  Company,  Cold  Spring,  Minnesota. 
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The  equations  above  are  based  on  the  assumptions  that 
(1)  the  material  is  linear  elastic,  (2)  lateral  inertia  forces 
are  negligible,  and  (3)  all  cross  sections  of  the  bar  remain 
plane  during  deformation. 

Tensile  strength  measurements  using  the  wave  propagation 
technique  depend  on  the  phenomenon  of  stress-wave  reflection  from 
a  free  boundary.  The  compressive  pulse  introduced  at  one  end  of 
the  bar  is  reflected  as  a  mirror-image  tensile  pulse.  Since  the 
tensile  strength  of  rocks  is  much  lower  than  the  compressive 
strength,  the  bar  fails  in  tension  upon  reflection  of  the  pulse 
from  the  free  end. 

Marble  and  granite  bars  were  instrumented  with  strain 
gages  at  two  stations  located  4  in.  from  each  end  (Fig.  12). 

The  gages  were  connected  to  potentiometric  circuits  and  their 
signals  recorded  on  an  oscilloscope.  Preliminary  tests  were 
conducted  by  dropping  a  steel  ball  on  the  end  of  the  bar  from 
heights  of  19  in.  and  38  in.  Typical  records  obtained  from  marble 
specimens  are  shown  in  Fig.  13.  These  signals  show  appreciable 
attenuation  and  wave  scattering  due  to  internal  friction  losses 
in  the  material.  Propagation  velocities  computed  from  these 
records  are  relatively  low,  the  highest  velocity  measured  being 
100,000  in. /sec.  Similar  tests  on  granite  bars  gave  wave  velocity 
values  of  the  order  of  160,000  in. /sec,  which  are  lower  than 
expected . 

The  specimens  were  subsequently  loaded  at  one  end  with 
explosive  charges  (R.  Stresau  Laboratories,  N-l  detonators, 
approximately  50  mg  total  charge).  Figure  14  shows  strain  gage 
data  from  one  of  these  tests  in  marble.  Appreciable  scatter 
was  noticed  in  the  measured  wave  propagation  velocities.  The 
average  of  five  measurements  gave  a  value  of 
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Fig.  13  STRAIN  GAGE  SIGNALS  IN  MARBLE  BAR  LOADED  BY  FALLING  WEIGHT 
38  ?nD>  (*J  HEIGHT  °F  DR0P  19  1". ,  (b)  HEIGOT  OF  DROP 
Lower’Tr^rf  TJ3Ce  i8/1™*  Gage  4  i"-  from  Loaded  End" 

Lower  Trace  ls  Second  Gage  4  in.  from  First  Gage) 
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STRAINS  IN  MARBLE  BAR  LOADED  EXPLOSIVELY  AT  ONE  EN 
(I'PPER  TRACE  IS  FIRST  CAGE  4  in.  FROM  LOADED  END- 
LOWER  TRACE  IS  SECOND  CAGE  3.91  in.  FROM  FIRST  GAG 
AND  4  in.  FROM  UNLOADED  END) 


35 


This  was  based  on  the  time  of  travel  of  the  compressive  pulse 
between  the  two  instrumented  stations.  Upon  reflection  from  the  free 
unloaded  end  the  wave  velocity  decreased  by  approximately 
50  percent.  The  dynamic  modulus  computed  from  Eq ,  (18)  and 
a  measured  value  for  the  density  is 

E  =  (155, 000)2  (250  x  10'6)  =  6.0  x  106  psi 

The  reflected  tensile  pulse  invariably  fractured  the 
specimen  before  reaching  its  peak  value,  thus  this  pulse  has 
always  a  lower  amplitude  than  the  incident  compressive  pulse. 

Only  the  portion  of  the  reflected  tensile  pulse  which  fractures 
the  specimen  is  transmitted  beyond  the  point  of  fracture.  Thus, 
by  measuring  the  amplitude  of  the  tensile  pulse  transmitted 
beyond  this  point,  one  can  determine  the  dynamic  t  nsile  strength 
of  the  material.  The  average  of  eight  determinations  of  the 
dynamic  tensile  strength  of  marble  was 

ofc  =  2200  psi 

Figure  15  shows  results  from  similar  tests  on  granite 
bars.  The  average  of  four  tests  resulted  in  a  uniaxial  wave 
propagation  velocity  of 

c  =  210,000  in. /sec 
o 

Upon  reflection  from  the  free  unloaded  end  the  wave  propagation 
velocity  decreased  by  approximately  15  percent.  The  dynamic 
modulus  computed  from  this  velocity  and  the  measured  density 
of  the  material  is 

E  -  (210, 000)2  (246  x  10"6)  =  10.8  x  106  psi 

As  in  the  case  of  marble,  the  reflected  tensile  pulse  fractured 
the  specimen  in  every  case.  The  tensile  strength  measurements 
showed  appreciable  scatter  due  to  the  high  inhomogeneity  of  the 
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material.  The  average  dynamic  tensile  strength  was 

at  =  5500  psi 

5 . 2  Wave  Propagation  in  Marble  Plates 

5.2.1  Photoelastic  Coating  Experiments 

The  specimens  used  in  this  study  were  10  in.  x  12  in. 
x  1/2  in.  Vermont  marble  plates.  These  plates  were  coated  with 
0.080  in.  thick  photoelastic  coatings  and  loaded  on  one  edge  with 
explosive  charges  between  100  mg  and  400  mg  (R.  Stresau  Labora¬ 
tories,  BM-9  and  BM-2  detonators). 

Photoelastic  fringe  patterns  on  the  opaque  rock 
specimens  were  recorded  with  a  Beckman  and  Whitley  camera  type 
189,  operating  at  rates  between  250,000  and  500,000  frames  per 
second  (Fig.  16).  This  camera  is  based  on  the  Miller  motion- 
compensation  technique  which  eliminates  the  relative  motion  between 
the  light  rays  reflected  from  a  rotating  mirrjr  and  the  stationary 
film.  The  concept  consists  of  a  stationary  objective  lens,  a  rotat 
ing  mirror,  and  one  pair  of  stationary  refocusing  lenses  for  each 
frame.  These  refocusing  lenses  are  located  on  arcs  of  a  circle 
whose  common  center  is  the  center  of  rotation  of  the  mirror.  A 
strip  of  stationary  film  is  arranged  along  the  arc  of  another 
circle,  also  with  the  mirror  as  center.  The  objective  lens 
focuses  a  stationary  primary  image  of  say  a  fixed  object  onto  or 
near  the  surface  of  the  rotating  mirror.  When  the  reflected  beam 
in  its  rotation  passes  over  a  typical  refocusing  lens  pair,  a 
stationary  image  is  formed  on  the  fixed  film.  Extremely  fast 
shuttering  is  effected  as  the  rotating  mirror  scans  the  light 
across  successive  pairs  of  refocusing  lenses,  thus  producing 
multiple  frames.  The  framing  rate  of  cameras  of  this  concept, 
as  with  that  of  a  smear  design,  is  directly  proportional  to  the 
rotational  speed  of  the  mirror.  This  camera  is  capable  of  record¬ 
ing  twenty-five  35mm  frames  at  up  to  2.5  million  frames  per  second. 
Numerous  light  sources,  such  as  ruby  lasers,  argon  flash  bombs,  etc 
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Fig.  16  BECKMAN  AND  WHITLEY  MODEL  189  CAMERA 


39 


are  available  with  this  camera.  Complete  circuitry  is  provided 

for  synchronizing  camera,  illumination  and  event.  The  model  at 

IITRI  has  both  framing  and  streak  (or  smear)  capability.  The 

conventional  framing  capability  of  such  models  has  been  used  by 
23 

Flynn  to  record  conventional  dynamic  photoelastic  (transmission) 
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fringe  photographs.  Cole  and  associates  have  employed  this 
type  of  camera  on  birefringent  coatings  to  record  (both  in  color 
and  black  and  white)  wave  propagation  in  aluminum  structures, 
while  deGraf^  used  it  to  study  fracture  in  steel.  Bieniawski^ 
has  used  a  similar  type  of  camera  to  study  photoelastically 
fracture  and  crack  propagation  in  rocks.  Rowlands,  Taylor  and 
Daniel  have  combined  pulsed-laser  illumination  with  the 

smear  capability  of  the  Beckman  and  Whitley  type  189  to  record 
photoelastically  stress  waves  and  fracture  in  polymeric  materials. 

The  most  critical  problem  in  application  of  this  camera 
was  the  light  source.  The  applicability  of  the  ruby  laser  was 
investigated,  but  the  pulsed  ruby  laser  available  was  not  suit¬ 
able  for  the  application.  Two  other  types  of  light  source  were 
used,  the  Argon  bomb  and  an  electronic  flash.  The  Argon  bomb 
used  was  a  balloon  approximately  12  in.  in  diameter  and  36  in. 
long  filled  with  Argon  gas.  This  gas  is  ignited  with  an  RP-1 
detonator  (Reynolds  Industries,  Inc.)  and  a  2  in.  x  2  in.  x  1/4  in. 
piece  of  Detasheet  (E.  I.  DuPont  and  Company)  detonated  with  a 
2000  v  pulse  from  the  Beckman  and  Whitley  camera.  A  Kociak  gela- 
filter  No,  55  was  used  in  front  of  the  objective  lens  of  the 
camera  to  isolate  a  band  of  light  around  the  5200  Angstrom  wave¬ 
length  . 

Satisfactory  isochromatic  fringe  patterns  were  obtained 
at  1,000,000  frames  per  second  when  one  Argon  bomb  was  used 
behind  a  transparent  Plexiglas  plate  loaded  with  a  detonator. 
However,  the  light  from  the  Argon  bomb  was  insufficient  when 
reflected  from  the  photoelastic  coating  applied  to  the  rock 
specimen.  To  avoid  complete  destruction  of  the  specimen  the 
Argon  bomb  was  directed  at  a  mirror  which  reflected  the  light 
to  the  specimen. 
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In  order  to  utilize  more  of  the  light  it  was  attempted 
to  record  photoelastic  coating  fringe  patterns  in  color  by  elimi¬ 
nating  the  monochromatic  filter  and  by  using  high-speed  color 
film  (High-Speed  Ektachrome  film  with  160  ASA  rating,  force- 
developed  to  650  ASA  and  1250  ASA)  ,  Although  this  approach  pro¬ 
duced  visible  fringe  patterns,  it  is  not  the  most  preferred  one 
because  of  the  graininess  of  the  film  and  the  less  than  faithful 
reproduction  of  colors. 

Another  approach  for  utilizing  more  of  the  available 
light  was  considered  by  improving  the  reflectivity  of  the  interface 
between  the  photoelastic  coating  and  the  rock  specimen.  In  the 
usual  applications,  the  photoelastic  coatings  are  bonded  to  the 
specimen  with  an  aluminum  filled  cement  which  provides  a  diffused 
reflective  surface  of  low  reflectivity.  To  improve  the  amount  of 
reflected  light,  transparent  photoelastic  coatings  were  vacuum- 
deposited  with  aluminum  to  provide  a  mirror  surface  and  were  then 
bonded  to  the  rock  specimen  with  clear  cement.  Several  tests  were 
conducted  with  such  specimens  using  the  light  of  two  Argon  bombs 
reflected  by  a  mirror,  but  the  results  were  unsatisfactory. 

A  new  light  source  was  prepared  for  use  with  the  Beckman 
and  Whitley  camera.  It  is  an  electronic  flash  unit  employing  a 
Xenon  flash  tube.  It  is  similar  to  the  Beckman  and  Whitley  Model 
357  un^-t  incorporating  some  modifications  proposed  by 

Flynn.  It  was  found  after  several  tests  that  the  light  produced 

by  one  flash  tube  was  comparable  to  that  of  one  Argon  bomb  and, 
therefore,  inadequate.  The  possibility  existed  of  adding  more  flash 
tubes  to  the  system  but  was  not  pursued  for  financial  reasons. 

The  effort  was  finally  directed  toward  multiple  Argon 
bombs  aimed  directly  at  the  specimen.  Photoelastic  coatings  with 
the  usual  aluminum  filled  cement  were  used.  A  variety  of  films 
and  developers  with  varying  development  times  were  used.  Films 
used  were  4X  developed  in  Acufine  or  D76  developer,  Tri-X  Pan 
developed  in  FG-7  developer,  Kodak  2485  film  developed  in  developer 
No.  857,  and  Double-X  developed  in  Dektol„ 
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figure  17  shows  one  experimental  setup  used  employing  two 
Argon  bombs  located  at  a  distance  of  39  in,  from  the  specimen 
off  the  centerline,  Sixteen  frames  of  a  twenty-five  frame  record 
of  isochromatic  fringe  patterns  in  the  photoelastic  coating  on  a 
marble  specimen  are  shown  in  Fig.  18  The  specimen  was  loaded 
on  the  edge  with  a  BM-9  detonator  (R,  Stresau  Laboratory;  main 
charge  0.111  gm  PETN ,  priming  charge  0.110  gm  Lead  Azide).  The 
camera  was  operated  at  a  rate  of  500,000  frames  per  second.  The 
film  used  was  4X  developed  in  Acufine  developer.  Figure  19  shows 
similar  fringe  patterns  obtained  with  a  marble  specimen  In 
this  case  four  Ar£on  bombs  were  used  as  a  light  source.  Tri-X 
Pan  film  developed  in  FG-7  developer  was  used. 

Better  results  were  obtained  using  Double-X  film  developed 
in  Dektol  developer  for  4-1/2  minutes.  Figures  20  and  21  show 
isochromatic  fringe  patterns  for  two  marble  specimens  loaded  with 
a  BM-9  detonator.  A  similar  marble  specimen  was  loaded  with  a 
BM-2  detonator  (R.  Stresau  Laboratory;  main  charge  0,331  gm  PETN, 
priming  charge  0,110  gm  Lead  Azide).  The  camera  was  operated  at 
a  rate  of  250,000  frames  per  second.  The  purpose  of  this  test 
was  to  span  a  longer,  up  to  100  \l  sec,  framing  time  and  observe 
a  later  stage  of  the  wave  propagation  phenomenon.  Channel  No.  1 
of  the  camera  was  used  to  trigger  the  detonator  and  two  Argon 
bombs  10  [i  sec  before  frame  No.  1,  Channel  No.  2  was  delayed  by 
59  p.  sec  with  respect  to  Channel  No0  1  and  was  used  to  detonate 
two  more  Argon  bombs.  Figure  22  shows  frames  16  through  23  of 
the  record  obtained. 

The  photoelastic  records  obtained  were  analyzed  to  de¬ 
termine  wave  propagation  and  attenuation  characteristics.  The 
propagation  velocities  of  various  fringe  orders  were  measured 
and  plotted  versus  fringe  order  and,  by  extrapolation  to  the 
zero  order  fringe,  the  wavefront  velocity  was  obtained.  Con¬ 
siderate  scatter  was  observed  in  the  results  from  many  specimens. 
Values  for  the  measured  dilatational  velocity  ranged  between 
148,000  and  172,000  in, /sec.  Typical  plots  of  propagation 
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Fig.  18  ISOCHROMATIC  FRINGE  PATTERNS  IN  PHOTOELASTIC  COATING  ON  MARBLE  SPECIMEN  LOADED 
EXPLOSIVELY  ON  THE  EDGE.  (CAMERA  SPEED:  500,000  FRAMES  PER  SECOND;  TIME  RA  GE 
14-44  M-  SEC;  SPECIMEN  100372-1) 
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velocity  versus  fringe  order  are  shown  in  Fig,  23  for  four 
specimens . 


In  addition  to  the  experiments  above,  marble  plates  were 
instrumented  with  strain  gages  and  loaded  on  the  edge  with 
explosive  charges  (BM-9  detonators) .  From  strain  gage  signals 
recorded  on  the  oscilloscope  a  plate  velocity  of  approximately 
160,000  in- /sec  was  determined. 

The  dilat&tional  velocity  in  a  plate  under  plane  stress 
conditions  is  given  by 


pd-vz) 


(20) 


Using  the  measured  wave  velocity  of  155,000  in- /sec  in  marble 
bars  and  a  value  for  v  =  0.25,  the  plate  velocity  is  computed 
as  c^  =  161,000  in. /sec. 

The  shear  wave  propagation  velocity  was  measured  from  the 
fringe  patterns  of  Fig.  22.  Because  of  the  finite  width  of  the 
plate  the  incident  pulse  is  reflected  from  the  free  sides  of  the 
specimen.  The  most  dominant  reflected  wave  is  the  shear  wave 
and  is  manifested  in  the  form  of  nearly  straight  parallel  fringes 
coming  in  from  the  sides  of  the  specimen.  The  velocity  of  these 
fringes  in  the  direction  normal  to  them  was  measured  and  was 
found  to  remain  nearly  constant  with  fringe  order.  The  measured 
shear  wave  propagation  velocity  was 

c  =  100,000  in. /sec 

The  shear  wave  velocity  was  also  computed  from  the  measured 
dilatational  velocity  and  Poisson's  ratio  using  relation  (17): 


c 


s 


160,000 


1-0-25 
- 2 


98,000  in. /sec 
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PROPAGATION  VELOCITY,  V,  ( in/sec) 


An  attempt  was  made  to  identify  a  Rayleigh  surface  wave 
along  the  top  edge  of  the  specimen.  An  area  of  birefringence 
concentration  on  the  edge  can  be  discerned  in  the  fringe  patterns 
of  Figs.  20  and  21.  The  propagation  velocity  of  this  pattern 
was  measured  to  be  65,000  in. /sec „ 


The  fringe  patterns  of  Fig.  20  were  further  analyzed  to 
determine  propagation  and  attenuation  characteristics.  Figure  24 
shows  the  variation  with  time  and  distance  of  pulse  shape  and 
amplitude  in  terms  of  fringe  order.  The  pulse  has  an  undulating 
shape  and  the  primary  leading  compressive  pulse  shows  appreciable 
attenuation  but  without  much  noticeable  dispersion.  To  separate 
the  effects  of  geometric  and  material  attenuation,  an  expression 
of  the  form 


n  e“ky 
n  =  J2 - 

max  Vy 


(21) 


for  the  maximum  fringe  order  was  sought,  where  y  is  the  distance 

from  the  loading  source  and  k,  the  coefficient  of  attenuation. 
-1/2 

The  term  y  accounts  for  geometric  attenuation  in  a  two- 

dimensional  space  (plate).  The  attenuation  coefficient  and  the 
constant  nQ  were  obtained  by  plotting  y1//2  versus  y  on  a 

semilog  scale  (Fig,  25).  Thus,  for  the  dilatational  wave  in 
Vermont  marble,  we  have 


n 


max 


4.78e~0,101y 

/T 


The  loss  tangent  of  the  material  for  a  sinusoidal  pulse  of 
wavelength  A  is  related  to  the  attenuation  coefficient  by 


tan  6 


(22) 


For  the  dilatational  wave  in  marble  the  loss  tangent  is  computed 
as 


tan  6 


^  0.101  x  3 

7T 


0.096 
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Fig.  25  SEMILOGARITHMIC  PLOT  FOR  SEPARATING  MATERIAL 
AND  GEOMETRIC  ATTENUATION 
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5.2.2.  Moir£  Experiments 


Marble  plates  similar  to  those  used  in  the  photoelastic 
coating  experiments  were  coated  with  1000  line-per-inch  rulings 
for  application  of  moiri  techniques.  The  specimen  rulings  were 
reproduced  on  0.004  in.  thick  Kodaline  reproduction  film  with 
estar  base,  and  the  film  was  bonded  onto  the  specimen  with  the 
emulsion  side  outwards  using  an  aluminum-filled  cement  (Photo- 
lastic  PC-8  wit n  PCH-8  hardener) .  The  purpose  of  the  aluminum 
filling  was  to  provide  a  somewhat  reflective  interface  for 
easier  photography  of  fringe  patterns.  The  specimen  rulings  were 
oriented  parallel  to  the  loaded  edge  of  the  specimen  to  yield 
full-field  representation  of  the  displacements  normal  to  this 
edge.  A  film  master  with  the  same  rulings  was  placed  in  contact 
and  aligned  with  the  specimen  ruling  by  means  of  a  film  of  oil. 

The  ruling  density  used,  representing  a  practical  upper 
limit  for  the  current  application,  is  not  very  sensitive  with  high 
modulus  materials  such  as  rock.  To  obtain  an  adequate  analyzable 
response  a  higher  explosive  loading  was  used.  This  consisted  of 
an  RP-1  detonator  (Reynolds  Industries,  Inc.,  250  mg  PETN  initiat¬ 
ing  explosive,  385  mg  Tetryl  high  density  explosive  and  4.1  gm 
RDX  booster  explosive),  and  two  3/4  in.  diameter  tetryl  pellets 
placed  under  the  first  detonator.  A  small  3/16  in.  thick  steel 
plate  was  placed  between  the  pellets  and  the  specimen.  Illumi¬ 
nation  was  provided  with  two  Argon  bombs  placed  at  distances 
between  25  in,  and  32  in.  from  the  specimens.  Moir£  fringe 
patterns  during  wave  propagation  were  photographed  with  the 
Beckman  and  Whitley  camera  at  rates  between  500,000  and 
1,000,000  frames  per  second.  Doule-X  film  developed  in  Dektol 
or  D-76  developer  was  used. 

Figure  26  shows  moir£  fringe  patterns  in  a  marble  specimen 
loaded  with  an  RP-1  detonator  and  two  tetryl  pellets  through  a 
3/16  in.  steel  plate.  The  two  Argon  bombs  were  placed  at  a 
distance  of  32  in.  from  the  specimen.  The  camera  was  operated 
at  a  rate  of  1,004,500  frames  per  second.  Figure  27  shows 
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similar  moir£  patterns  for  a  similarly  loaded  marble  specimen. 
The  Argon  bombs  were  placed  at  a  distance  of  25  in.  from  the 
specimen  and  the  camera  was  operated  at  a  rate  of  508,000  frames 
per  second. 


The  moir£  fringes  obtained  were  interpreted  and  analyzed 
as  loci  of  points  having  the  same  component  of  displacement 
normal  to  the  ruling  lines  of  the  undistorted  reference  (master) 
array.  Displacement  profiles  as  a  function  of  distance  from 
the  source  for  the  two  specimens  above  with  time  (or  frame 
number)  as  a  parameter  are  shown  in  Figs.  2f  and  29.  The  atten¬ 
uation  and  dispersion  of  the  displacement  pulse  is  evident. 

The  pulse  even  shows  sign  reversal  at  later  times  (Frames  20-24 
in  Fig.  29) . 

The  vertical  displacements  above  are  sufficient  for  the 
complete  determination  of  stresses  along  the  vertical  centerline 
for  the  case  of  a  dilatational  wave.  The  strains  along  the 
y-axis  are  related  to  the  displacements  as  follows: 


C  =  Jx 
-y  °y 


(23) 


from  which  the  stresses  or  and  o  can  be  determined.  The 

x  y 

vertical  strains  determined  from  the  displacements  of  Fig.  29 
are  plotted  in  Fig.  30  as  a  function  of  distance  from  the  source 
with  time  as  a  parameter.  It  is  seen  that  the  vertical  strain 
becomes  tensile  at  the  trailing  end  of  the  pulse  between  frames 
14  and  16.  Since  for  the  most  part  of  the  phenomenon  the  trans¬ 
verse  strains  are  positive  it  follows  that  tensile  stresses  are 
generated  at  the  trailing  end  of  the  pulse.  Given  the  low 

tensile  strength  of  rock,  these  tensile  stresses  are  very  likely 
to  produce  fracture. 
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Fig.  28  VERTICAL  DISPLACEMENT  ALONG  Y-AXIS  AS  A  FUNCTION  OF  DEPTH 

A  PARAMETER  (FRAME  NUMBERS  MARKED;  CAMERA 
SPEED:  1,004,500  FRAMES/SEC;  MARBLE  SPECIMEN  41673-2) 
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EXPLOSIVE 


5 . 3  Wave  Propagation  in  Granite  Plates 


The  specimens  used  were  10  in.  x  12  in.  x  1/2  in.  Cold 
Spring  (Minnesota)  granite  plates.  They  were  coated  with  0.080  in. 
thick  photoelastic  coatings.  They  were  loaded  explosively  on 
one  edge  and  the  resulting  isochromatic  fringe  patterns  photo¬ 
graphed  with  the  Beckman  and  Whitley  camera  operating  at  a  rate 
of  500,000  frames  per  second. 

Figure  31  shows  fringe  patterns  obtained  from  one  granite 
specimen  loaded  with  a  BM-9  detonator.  Two  Argon  bombs  at  a 
distance  of  39  in.  from  the  specimen  were  used  as  a  light  source. 
4X  film  was  used  and  developed  in  Acufine  developer.  Figure  32 
shows  similar  fringe  patterns  obtained  from  a  similar  granite 
specimen  loaded  with  a  BM-2  detonator  on  the  edge.  Three  Argon 
bombs  placed  at  a  distance  of  36  in.  from  the  specimen  were  used 
as  the  light  source.  Double-X  film  processed  in  Dektol  developer 
was  used.  The  fringe  patterns  -in  both  cases  are  not  well  defined 
because  of  the  pronounced  inhomogeneity  of  granite.  The  material 
seems  to  consist  of  large  grains  of  varying  mechanical  and 
thermal  properties.  Temperature  changes  from  the  temperature  at 
which  the  coating  was  bonded  induce  nonuniforn;  fringe  patterns 
which  are  superimposed  on  the  patterns  induced  by  loading. 

The  wave  propagation  of  various  fringe  orders  was 
measured  and  plotted  versus  fringe  order  in  Fig.  33.  The  dilata- 
tional  wave  propagation  velocity  obtained  by  extrapolation  to 
the  zero-order  fringe  was  172,000  and  192,000  in. /sec  for  the 
two  specimens.  As  in  the  case  of  marble  appreciable  attenuation 
is  evident. 
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6-0  FRACTURE  PROPAGATION  IN  ROCK  MEDIA 

It  was  shown  earlier  that  the  initially  compressive  pulse 
undergoes  sign  changes  with  tensile  stresses  developing  at  the 
trailing  end  at  later  times.  If  the  initial  pulse  is  sufficiently 
intense  the  geometric  and  material  attenuation  will  not  reduce 
these  tensile  stresses  to  a  level  below  the  tensile  strength  of 
the  material.  This  phenomenon  results  in  a  fracture  zone  around 
the  point  source  clearly  evident  as  the  nearly  semicircular 
region  of  dense  fringes  in  Figs.  20-22.  The  propagation  velocity 
of  the  boundary  of  this  zone  was  measured  in  the  marble  specimens 
of  Figs.  20  and  21  as  30,000  and  35,000  in. /sec.  A  similar 
fracture  zone  develops  in  granite  specimens.  These  fractures 
are  illustrated  in  Fig.  34  for  marble  and  granite.  In  the  case 
of  the  granite  specimen  additional  fractures  near  the  upper 
corners  were  caused  by  tensile  pulses  reflected  from  the  sides 
of  the  specimen. 

Several  attempts  were  made  to  induce  crack  propagation 
in  the  rock  plates.  A  notch  was  introduced  at  the  bottom  of  the 
specimen  opposite  the  loaded  edge  with  the  expectation  that  the 
tensile  pulse  reflected  from  the  bottom  edge  would  initiate  a 
crack.  This,  however,  did  not  occur,  probably  because  of  the 
high  attenuation  of  the  wave.  A  notch  was  then  introduced  at 
the  upper  edge  and  the  specimens  loaded  with  explosive  packed 
inside  this  notch.  Figure  35  shows  isochromatic  fringe  patterns 
obtained  from  a  marble  specimen  loaded  explosively  through  a 
notch.  An  RP-1  detonator  and  2.23  gm  of  Deta  sheet  were  used. 

Four  Argon  bombs  placed  at  a  distance  of  36  in.  from  the  specimen 

were  used  as  a  light  source.  The  film  was  Double-X  developed  in 
Dektol . 

The  propagation  of  the  fracture  zone  identified  as  the 
area  of  dense  fringe  pattern  was  measured  (Fig.  36).  The  velocity 

of  propagation  of  the  boundary  of  this  zone  is  not  uniform  but  a 
value  of 

vfracture  =  50,000  in. /sec 
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Fig.  36  PROPAGATION  OF  FRACTURE  ZONE  BOUNDARY  IN  MARBLE 
PIRATE  LOADED  EXPLOSIVELY  THROUGH  A  NOTCH 
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was  measured  near  the  central  part  of  the  curve. 

Additional  marble  and  granite  plates  were  loaded  ex¬ 
plosively  with  an  RP-1  detonator  and  a  tetryl  pellet  placed  over 
a  steel  wedge  inserted  in  a  notch  on  one  edge.  Four  Argon  bombs 
and  Double-X  film  were  used.  The  camera  was  operated  at  a  rate 
of  500,000  frames  per  second.  Figures  37  and  38  show  isochromatic 
fringe  patterns  in  the  photoelastic  coating  around  propagating 
cracks  in  the  marble  and  granite  plates.  The  crack  propagation 
velocities  measured  were  38,000  in. /sec  and  36,000  in„/sec  for 
the  marble  and  granite  specimens,  respectively. 
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Fig.  37  CRACK  PROPAGATION  IN  MARBLE  PLATE  LOADED  EXPLOSIVELY  THROUGH  METAL  WEDGE 
IN  NOTCH  (FRAMES  SHOWN  AR1  AT  4  |i  SEC  INTERVALS;  SPECIMEN  31573-3) 


7.0  SUMMARY  AND  DISCUSSION 

Experimental  stress  analysis  techniques,  such  as  dynamic 
photoelasticity,  moir£  and  photoelastic  coatings  were  used  to 
study  wave  and  fracture  propagation  in  rock  and  rocklike  media. 
Materials  used  in  t  is  investigation  were  transparent  rocklike 
birefringent  materials  such  as  Columbia  Resin  CR-39,  Plexiglas 
and  Homalite-100  and  two  rock  materials,  namely,  Vermont  Marble 
and  Cold  Spring  Granite. 

The  transparent  birefringent  specimens  were  loaded  with 
explosive  charges  at  an  interior  point  and  the  isochromatic 
fringe  patterns  and  crack  propagation  were  recorded  with  a 
Cranz-Schardin  multiple  spark  camera  operating  at  a  rate  of 
200,000  frames  per  second.  The  photographic  records  obtained 
were  analyzed  to  determine  dilatational  wave,  shear  wave  and 
crack  propagation  velocities.  Some  of  the  results  are  sum¬ 
marized  in  Table  2.  In  the  cast  of  the  brittle  CR-39  material, 
it  was  possible  to  measure  crack  propagation  velocities  up  to 
half  the  measured  shear  wave  velocity,  which  represents  the 
theoretical  limit  for  fracture  propagation.  No  fracture  propa¬ 
gation  was  observed  in  Plexiglas  and  limited  fracture  was 
observed  in  Homalite-100.  The  crack  propagation  velocity 
measured  in  the  latter  was  only  0.273  times  the  shear  wave 
velocity,  far  below  the  theoretically  predicted  limiting 
velocity. 

Marble  and  granite  plates  were  loaded  explosively  on 
one  edge  and  isochromatic  fringe  patterns  on  bonded  photoelastic 
coatings  and  moiri  fringe  patterns  were  recorded  with  a  Beckman 
and  Whitley  camera  operating  at  rates  from  250,000  to  1,000,000 
frames  per  second. 
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Dilatational,  shear  and  Rayleigh  wave  velocities  were 
measured  in  marble.  Some  of  the  results  are  tabulated  in  Table 
2.  The  measured  shear  wave  velocity  was  almost  identical  to 
the  one  computed  from  the  measured  dilatational  wave  velocity 
and  Poisson's  ratio  of  the  material.  The  characteristics  of 
the  propagating  pulse  were  studied.  The  primary  leading  part 
of  the  pulse  is  compressive  and  shows  appreciable  attenuation 
but  not  much  noticeable  dispersion.  The  attenuation  observed 
was  separated  into  its  geometric  and  material  components.  The 
coefficient  of  material  attenuation  in  marble  was  found  to  be 
k  ■  0.101  which  corresponds  to  a  loss  tangent  of  tan  6  =  0.096. 

Transient  displacement  and  strain  distributions  were 
determined  from  the  moir£  fringe  patterns.  It  was  shown  that 
the  trailing  end  of  the  propagating  pulse  becomes  tensile,  thus 
likely  to  induce  widespread  tensile  fracture.  The  dilatational 
wave  propagation  velocity  was  also  measured  in  granite,  which 
also  exhibited  appreciable  attenuation. 

To  study  fracture  propagation  marble  and  granite  plates 
were  loaded  explosively  through  a  notch  on  one  edge.  The  crack 
propagation  velocities  measured  were  relatively  low,  below  the 
limiting  theoretical  values.  Crack  propagation  velocities  of 
38,000  in. /sec  and  36,000  in. /sec  were  determined  for  marble 
and  granite,  respectively.  The  corresponding  theoretical 
terminal  velocities  are  50,000  and  56,000  in. /sec.  In  the  case 
of  one  marble  specimen  loaded  with  a  heavy  explosive  inside  the 
notch,  the  propagation  velocity  of  the  fracture  ’jot-c  was 
measured  to  be  50,000  in. /sec.  It  seems  that  terminal  fracture 
velocities  in  rock  can  be  achieved  if  the  pulse  is  sufficiently 
intense.  The  high  material  and  geometric  attenuation  tends  to 
weaken  the  propagating  pulse  and  reduce  the  fracture  propagation. 
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The  experimental  investigation  conducted  and  described 
above  was  of  a  fundamental  nature,  aimed  at  providing  a  better 
understanding  of  the  rapid  excavation  process  by  explosives. 

The  primary  input  of  an  explosive  in  a  rock  mass  is  a  compressive 
pulse  traveling  with  the  dilatational  wave  velocity.  The  geo¬ 
metric  attenuation  of  such  a  pulse  is  very  rapid  in  three- 
dimensional  space,  the  pulse  amplitude  varying  inversely  with  the 
distance  from  the  explosive  source.  This  fact,  coupled  with  the 
high  compressive  strength  of  rock,  tends  to  limit  compressive 
fracturing  in  the  immediate  vicinity  of  the  explosive  source. 
However,  high  speed  explosives  tend  to  produce  short  compressive 
pulses  followed  by  tensile  trailing  pulses,  as  evidenced  in  the 
model  studies  described  in  this  report.  The  resulting  tensile 
stresses,  even  though  they  are  of  relatively  low  magnitude,  tend 
to  produce  more  extensive  tensile  fracture  after  passage  of  the 
pulse,  because  of  the  much  lower  tensile  strength  of  rock. 

Tensile  rock  fracture  also  results  when  the  primary  com¬ 
pressive  pulse  interacts  with  preexisting  cavities  or  flaws  in 
the  rock  mass  and  when  it  is  reflected  as  a  tensile  pulse  from 
a  free  surface.  It  is  conceivable  that  the  efficiency  of  explosive 
excavation  could  be  enhanced  by  preinducing  such  cavities  or 
flaws  V7ith  prior  smaller  explosives  properly  spaced  and  timed. 

Multiple  explosives  suitably  timed  and  spaced  can  also 
be  used  to  produce  pulses  interfering  constructively  in  specified 
regions  of  the  rock  mass,  thus  producing  the  maximum  fracture 
for  a  given  total  amount  of  explosive. 

Multiple-explosive  drilling  in  rock  is  a  common  practice. 
Systems  have  been  used  whereby  explosive  charges  are  delivered 
to  the  bottom  of  a  well  at  a  steady  rate  and  detonated  sequenci- 
ally  while  at  the  same  time  flushing  up  the  fractured  rock  chips. 

A  variation  of  this  general  approach  has  been  developed  recently 
whereby  the  explosive  charge  consists  of  a  series  of  smaller 
conical ^charges  aimed  in  different  directions  away  from  the 
center.  This  device  thus  sends  out  a  series  of  small  jets  into 


1 1 T  RESEARCH  INSTITUTE 

75 


the  rock  in  a  conical  pattern.  It  is  claimed  that  this  system 
can  drill  more  rock  per  unit  of  explosive  than  previously  used 
systems . 

Further  developments  of  this  nature  can  definitely 
benefit  from  fundamental  laboratory  experiments  such  as  described 
here.  Knowledge  of  the  dynamic  rock  properties  and  the  pulse 
propagation  and  attenuation  characteristics  of  the  explosive 
would  help  in  properly  designing  a  series  of  multiple  explosions 
for  maximum  efficiency.  Further  experimental  laboratory  work 
with  rock  specimens  is  recommended.  High-speed  photographic 
techniques  can  be  improved  further  to  record  photoelastic  and 
moir£  fringe  patterns  for  a  more  detailed  analysis  of  experi¬ 
mental  results.  Means  for  inducing  crack  propagation  at  near 
the  theoretical  limiting  velocity  must  be  explored.  The  effects 
of  multiple  explosive  loadings  must  be  studied  to  determine  the 
optimuu  spacing  and  timing  of  these  explosives. 
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